APPLIED PHYSICS LETTERS xx (2011) 



o 

(N 
C 



Hybrid spintronics and straintronics: A magnetic technology for ultra low 
energy computing and signal processing 

Kuntal Roy/'EIS upriyo Bandyopadhyay, 1 and Jayasimha Atulasimha 2 

^Department of Electrical and Computer Engineering, Virginia Commonwealth University, Richmond, VA 23284, 
USA 

2 ' Department of Mechanical and Nuclear Engineering, Virginia Commonwealth University, Richmond, VA 23284, 
USA 

(Dated: 13 January 2011) 

The authors show that the magnetization of a magnetostrictive/piezoelectric multiferroic single-domain shape- 
anisotropic nanomagnet can be switched with very small voltages that generate strain in the magnetostrictive 
layer. This can be the basis of ultralow power computing and signal processing. With appropriate material 
choice, the energy dissipated per switching event can be reduced to ^45 kT at room temperature for a 
switching delay of ^100 ns and ~70 kT for a switching delay of ~10 ns, if the energy barrier separating the 
two stable magnetization directions is ^32 kT. Such devices can be powered by harvesting energy exclusively 
from the environment without the need for a battery. 
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The primary obstacle to continued downscaling of dig- 
ital electronic devices in accordance with Moore's law is 
the excessive energy dissipation that takes place in the 
device during switching of bits. Every charge-based de- 
vice (e.g. MOSFET) has a fundamental shortcoming in 
this regard. They are switched by injecting or extracting 
an amount of charge AQ from the device's active region 
with a potential gradient AV, leading to an inevitable 
energy dissipation of AQ x AV. Spin based devices, on 
the other hand, are switched by flipping spins without 
moving any charge in space (AQ = 0) and causing a 
current flow. Although some energy is still dissipated in 
flipping spins, it can be considerably less than the energy 
AQ x AV associated with current flow. This gives "spin" 
an advantage over "charge" as a state variable to encode 
digital bits. 

Recently, it has been shown that the minimum energy 
dissipated to switch a charge-based device like a transis- 
tor at a temperature T is ~NkTln(l/p), where N is the 
number of information carriers (electrons or holes) in the 
device and p is the bit error probability- 1 -. On the other 
hand, the minimum energy dissipated to switch a single 
domain nanomagnet (which is a collection of M spins) 
can be only ~fcTZn(l/p) since the exchange interaction 
between the spins makes all of them rotate together in 
unison like a giant classical spin 1 ? 2 . This gives the mag- 
net an advantage over the transistor. 

Unfortunately, the magnet's advantage is lost if the 
method adopted to switch it is so inefficient that the en- 
ergy dissipated in the switching circuit far exceeds the 
energy dissipated in the magnet. Regrettably, this is of- 
ten the case. A magnet is usually flipped with either 
a magnetic field generated by a current 3 -, or a spin po- 
larized current exerting either a spin transfer torque^ or 
causing domain wall motion 5 -. The energy dissipated to 
switch a magnet with current-generated magnetic field 
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was reported in ref. [3] as 10 11 - 10 12 kT for a switch- 
ing delay of ~1 fis, which clearly makes it impractical. 
On the other hand, the energy dissipated to switch a 
two-dimensional nanomagnet in ^1 ns with spin transfer 
torque is estimated to be ~2xl0 8 kT if the energy bar- 
rier separating the two stable magnetization orientations 
is 32 kT (so that the equilibrium bit error probability p 
= e~ 32 )£. Therefore, both of these switching methods 
are very energy-inefficient since the switching circuit dis- 
sipates far more energy than the minimum ~kTln(l/p) 
needed. In fact, they are so inefficient that they might 
not even make the magnet superior to the transistor 
which can be switched in sub-ns while dissipating 10 7 

- 10 8 kT of energy in a circuit. Only domain wall mo- 
tion has been shown to be a relatively energy-efficient 
switching mechanism since there is at least one report 
of switching a nanomagnet in 2 ns while dissipating 10 4 

- 10 5 kT of energy- 8 -. Thus, there is a need to identify 
energy-efficient mechanisms for switching a magnet. This 
is the motivation for this work. 

Recently, we have shown that the magnetization of 
a strain-coupled piezoelectric/magnetostrictive multifer- 
roic nanomagnet can be switched by stressing the mag- 
netostrictive layer with a small voltage applied to the 
piezoelectric layer—. Such multiferroic systems have now 
become commonplac e 10 ' 11 ! 13 and there are proposals for 
using them in magnetic logic and memor y 6 ' 9 ' 14 ' 15 . In this 
method, an electrostatic potential applied to the piezo- 
electric layer of a multiferroic magnet generates in it a 
strain that is elastically transferred to the magnetostric- 
tive layer if the latter layer is considerably thinner. This 
stresses the magnetostrictive layer and causes its mag- 
netization to rotate. Such rotations have been demon- 
strated experimentally^. Consider an ellipsoidal multi- 
ferroic magnet with uniaxial shape anisotropy as shown 
in Fig. 1. The magnetostrictive layer is assumed to be 
10 nm thick and the piezoelectric layer is 40 nm thick, 
which ensures that most of the strain generated in the 
piezoelectric layer by an applied voltage is transferred to 
the magnetostrictive layer. We assume that the piezo- 
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FIG. 1. An elliptical multiferroic nariomagnet stressed with 
an applied voltage. 



electric layer is lead-zirconate-titanate (PZT) and the 
magnetostrictive layer is polycrystalline nickel or cobalt 
or Terfenol-D. For Terfenol-D, the major axis is assumed 
to be ^102 nm and the minor axis is ~98 nm. Because 
of shape anisotropy, the two magnetization orientations 
parallel to the easy axis (major axis of the ellipse, or 
the z-axis) are stable and store the binary bits and 1. 
The energy barrier between these two orientations (i.e. 
the shape anisotropy barrier) will be 0.8 eV or ^32 kT 
at room temperature for the Terfenol-D system, which 
makes the equilibrium bit error probability e -32 . Let us 
assume that the magnetization is initially oriented along 
the -z-axis. Our task is to switch the nanomagnet so 
that the final orientation is along the +z-axis. We do 
this by applying a stress along the easy axis (z-axis). 
Since the stress is generated with a voltage V applied to 
the piezoelectric layer, the energy dissipated during turn- 
on is (1/2)CV 2 while that dissipated during turn-off is 
(1/2)CV 2 , where C is the capacitance of the piezoelectric 
layer, and we have assumed that the voltage is applied 
abruptly or non-adiabatically. 

There is an additional dissipation Ed in the nanomag- 
net due to Gilbert damping 12 , which is of the order of 
the energy barrier modulation in the nanomagnet, i.e. 
by how much the energy barrier between the two stable 
magnetizations is lowered to switch between the states. 
The total energy dissipated in the switching process is 
therefore E to tai = CV 2 + Ed- Thus, in order to calculate 
Etotai as a function of switching delay, we have to calcu- 
late four quantities: (1) the stress needed to switch the 
magnetization within the given delay, (2) the voltage V 
needed to generate this stress, (3) the capacitance C of 
the multiferroic, and (4) Ed which is calculated by follow- 
ing the prescription of ref. [12] (see the supplementary 
material) . 

In order to find the stress a required to switch a magne- 
tostrictive nanomagnet in a given time delay r, we solve 
the Landau-Lifshitz-Gilbert (LLG) equation for a single 
domain magnetostrictive nanomagnet subjected to stress 
a. We then relate a to the strain e in the nanomagnet 
from Hooke's law (e = a/Y, where Y is the Young's 
modulus of the nanomagnet) and find the voltage V that 
generates that strain in the PZT layer based on the d^i 
coefficient of PZT and its thickness. Finally, we calculate 
the capacitance of the multiferroic system by treating it 
as a parallel-plate capacitor while taking the relative di- 
electric constant of PZT to be 1000. This allows us to 
find the energy dissipated in the switching circuit (CV ) 



as a function of the switching delay t. 

In the supplementary material accompanying this let- 
ter, we show that both stress and shape anisotropy act 
like a torque on the magnetization of the nanomagnet. 
This torque per unit volume of the nanomagnet is given 
by 



Te(*) 



x VE[6(t),<p(t)] 



(1) 



where E[0(t),<j)(t)] is the total energy of the nanomagnet 
at an instant of time t. It is the sum of shape anisotropy 
energy and stress anisotropy energy, both of which de- 
pend on the magnetization orientation at the given in- 
stant. We adopt the spherical coordinate system whereby 
the magnetization is along the radial direction so that its 
orientation at any instant is specified by the instanta- 
neous polar angle 0(t) and the azimuthal angle 4>(t). 

In the supplementary material, we show that we can 
write the torque as 

T E (*) = -{2B( ( f)(t))sin0(t)cos0(t)}e < f, 

-{B Oe (^(t))sin0(t)}eg 7 (2) 

where eg and are unit vectors in the 0- and in- 
directions, and 

MM)) = y M s Q [N xx cos 2 <l)(t) + N yv sin 2 4>(t) - N zz ] 

(3a) 
(3b) 
(3c) 

(3d) 



Bstress — (3/2)A s aO 

B(4>(t)) = B (<l>(t)) + B a tress 



BoeW)) 



Mo 



M*Sl(N xx - N vv )sin(2<j>(t)). 



Here M s is the saturation magnetization of the nanomag- 
net, is its volume, /j,q is the permeability of free space, 
A s is the magnetostrictive coefficient of the magnetostric- 
tive layer, and Npp is the demagnetization factor in the 
ft direction, which can be calculated from the shape and 
size of the nanomagnet [see the supplementary material] . 

The magnetization dynamics of the single-domain 
nanomagnet is described by the Landau-Lifshitz-Gilbert 
(LLG) equation as follows. 



dn m (t) 
dt 



a n 



dn m (t) 
dt 



7 
M v 



T B (t) (4) 



where n m (i) is the nomalized magnetization, a is the di- 
mensionless phenomenological Gilbert damping constant, 
7 = 2/iB/xo/fi.is the gyromagnetic ratio for electrons, and 
M v = n M s Q. 

From this equation, we can derive two coupled equa- 
tions that describe the 0- and 0-dynamics. The deriva- 
tion can be found in the supplementary material. The 
final result is: 

7 

' M v ' 

+2aB(4>(t))sin0(t)cos0(t)}, (5) 
(1 + a 2 ) <j>'(t) = -1- [aB 0e (<f>(t)) - 2B(<t>(t))cos0(t)] 

(sin0(t)^O). (6) 



(1 
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[B 0e (<t>(t))sin6(i) 
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FIG. 2. Energy dissipated in the switching circuit (CV 2 ) and 
the total energy dissipated (E to tai) as functions of delay for 
three different materials used as the magnetostrictive layer in 
the multiferroic nanomagnet. 



Clearly, the 9- and ^-motions are coupled and hence 
these equations have to be solved numerically. We as- 
sume that the initial orientation of the nanomagnet is 
close to the -z-axis (9 — 179°). It cannot be exactly 
along the -z-axis (9 — 180°) since then the torque acting 
on it will be zero (see Equation ([2])) and the magneti- 
zation will never rotate under any stress. Similarly, we 
cannot make the final state align exactly along the +z- 
axis (9 = 0°) since there too the torque vanishes. Hence 
we assume that the final state is 9 — 1°. Thus, both ini- 
tial and final states are 1° off from the easy axis. Thermal 
fluctuations can easily deflect the magnetization by 1°. 

We apply the voltage generating stress abruptly at 
time t = 0. This rotates the magnetization away from 
near the easy axis (9 — 179°) since the latter is no longer 
the minimum energy state. The new energy minimum is 
at 9 = 90°. We maintain the stress until 9 reaches 90° 
which places the magnetization approximately along the 
in-plane hard axis (y-axis). Then we reduce the voltage 
to zero abruptly. Subsequently, shape anisotropy takes 
over and the magnetization vector will rotate towards the 
easy axis since that becomes the minimum energy state. 
The question is which direction along the easy axis will 
the magnetization vector relax to. Is it the -z-axis at 
9 = 179° (wrong state) or the +z-axis at 9 = 1° (correct 
state)? That is determined by the sign of Bo e (<fi(t)) when 
9 reaches 90°. If cj> at that instant is less than 90°, then 
Bo e ((j)(t)) is positive which makes the time derivative of 
9 negative (see Equation ((SJ) so that 9 continues to de- 
crease and the magnetization reaches the correct state 
close to the +z-axis. The coupled 9- and ^-dynamics en- 
sures that this is the case as long as the stress exceeds 
a minimum value. Thus, successful switching requires a 
minimum stress. 

Once we have found the switching delay r for a given 
stress a by solving Equations (J5)) and ([6]), we can invert 
the relationship to find a versus r and hence the energy 
dissipated versus r. This is shown in Fig. 2 where we 
plot the energy dissipated in the switching circuit (CV 2 ), 
as well as the total energy dissipated (E tota i) versus delay 
for three materials chosen for the magnetostrictive layer. 
For Terfenol-D, the stress required to switch in 100 ns 



is 1.92 MPa and that required to switch in 10 ns is 2.7 
MPa. 

Note that for a stress of 1.92 MPa, the stress anisotropy 
energy B stress is 32.7 kT while for 2.7 MPa, it is 46.2 kT. 
As expected, they are larger than the shape anisotropy 
barrier of ~32 kT which had to be overcome by stress to 
switch. A larger excess energy is needed to switch faster. 
The energy dissipated and lost as heat in the switching 
circuit (CV 2 ) is only 12 kT for a delay of 100 ns and 23.7 
kT for a delay of 10 ns. The total energy dissipated is 45 
kT for a delay of 100 ns and 70 kT for a delay of 10 ns. 
Note that in order to increase the switching speed by a 
factor of 10, the dissipation needs to increase by a factor 
of 1.6. Therefore, dissipation increases sub-linearly with 
speed, which bodes well for energy efficiency. 

With a nanomagnet density of 10 10 cm -2 in a mem- 
ory or logic chip, the dissipated power density would have 
been only 2 mW/cm 2 to switch in 100 ns and 30 mW/cm 2 
to switch in 10 ns, if 10% of the magnets switch at any 
given time (10% activity level). Note that unlike tran- 
sistors, magnets have no leakage and no standby power 
dissipation, which is an important additional benefit. 

Such extremely low power and yet high density mag- 
netic logic and memory systems, composed of multifer- 
roic nanomagnets, can be powered by existing energy 
harvesting systems^— that harvest energy from the 
environment without the need for an external battery. 
These processors are uniquely suitable for implantable 
medical devices, e.g. those implanted in a patient's brain 
that monitor brain signals to warn of impending epilep- 
tic seizures. They can run on energy harvested from 
the patient's body motion. For such applications, 10- 
100 ns switching delay is adequate. Speed is not the 
primary concern, but energy dissipation is. These hy- 
brid spintronic/straintronic processors can be also in- 
corporated in "wrist-watch" computers powered by arm 
movement, buoy-mounted computers for tsunami moni- 
toring (or naval applications) that harvest energy from 
sea waves, or structural health monitoring systems for 
bridges and buildings that are powered solely by mechan- 
ical vibrations due to wind or passing traffic. 
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In this supplementary section, we first derive the equations describing the time evolution of the 
polar angle 9{t) and the azimuthal angle <j){t) of the magnetization vector. We do this starting from 
the Landau-Lifshitz-Gilbert (LLG) equation. 

SI Magnetization dynamics of a multiferroic nanomagnet: Solu- 
tion of the Landau-Lifshitz-Gilbert equation 

Consider an isolated nanomagnet of ellipsoidal shape lying in the y-z plane with its major axis 
aligned along the z-direction and minor axis along the y-direction. The dimension of the major 
axis is a and that of the minor axis is b, while the thickness is I. The volume of the nanomagnet 
is = (7r/4)aM. Let 6(t) be the angle subtended by the magnetization with the +z-axis at any 
instant of time t and 4>{t) be the angle subtended by the projection of the magnetization vector on 
the x-y plane with the +x axis. We call 6{t) the polar angle and 4>(t) the azimuthal angle. These 
are represented in Fig. 1 of the main paper. 

The total energy of the single-domain nanomagnet is the sum of the uniaxial shape anisotropy 
energy E$ha and the stress anisotropy energy E$ta'- 



E = Esha + Est a, 



(SI) 
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where 



Esha = (p /2)M*nN d , 



(S2) 



with M s being the saturation magnetization and Nd the demagnetization factor expressed as 



N d = N zz cos 2 9(t) + N yy sin 2 6(t) sin 2 4>(t) + N xx sin 2 6(t) cos 2 (p(t) (S3) 



Here N zz , N yy , and N xx are the components of Nd along the z-axis, y-axis, and x-axis, respectively. 
If I <C a, b, then N zz , N yy , and N xx are given by [£Q] 
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(S4a) 

(S4b) 
(S4c) 



which shows that N xx 3> N yy , A^ 22 . 

Note that in the absence of any stress, uniaxial shape anisotropy will favor lining up the mag- 
netization along the major axis (z-axis) [9 = 0, (f> = 90°] by minimizing Esha, which is why we will 
call the major axis the "easy axis" and the minor axis (y-axis) the "hard axis". We will assume 
that a force along the z-axis (easy axis) generates stress in the magnet. In that case, the stress 
anisotropy energy is given by 

E S ta = -(3/2)\ s aQcos 2 6(t), (S5) 

where (3/2)A s is the magnetostriction coefficient of the nanomagnet and a is the stress. Note that a 
positive \ s a product will favor alignment of the magnetization along the major axis (z-axis), while a 
negative A s cr product will favor alignment along the minor axis (y-axis), because that will minimize 
Est A- In our convention, a compressive stress is negative and tensile stress is positive. Therefore, 
in a material like Terfenol-D that has positive X s , a compressive stress will favor alignment along 
the minor axis, and tensile along the major axis. The situation will be opposite with nickel that 
has negative X s . 

At any instant of time, the total energy of the nanomagnet can be expressed as 



E(t) = E[9{t),<j)(t)] = B(4>{t))sin 2 9{t) + C 



(S6) 
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where 



Mo 

Bo(Ht)) = — M 2 tt [N xx cos 2 0(t) + N yy sin 2 <j){t) - N zz ] (S7a) 
B s tress = (3/2)X s an (S7b) 

B(0(t)) = B (^(t)) + Stress (S7C) 

Mo „ 

C = —M 2 ttN zz - (3/2)A s fffi. (S7d) 

Note that Bo((/>(t)) is always positive, but B stress can be negative or positive according to the sign 
of the X s a product. 

The magnetization M(t) of the magnet has a constant magnitude at any given temperature but 
a variable direction, so that we can represent it by the vector of unit norm n m (t) = M(t)/|M| = e r 
where e r is the unit vector in the radial direction in spherical coordinate system represented by 
(r,8,<t>). The other two unit vectors in the spherical coordinate system are denoted by e# and for 
6 and 4> rotations, respectively. Note that 

8E(t) 1 dE(t) 

vE(t ) = vE m , m = + m 

dE{t) 

— — = 2Bsin6(t)cos9(t) (S9) 
ou(t) 

= M 2 fl(N xx - N yy )sin{2(p{t))sin 2 e{t) = -B Qe {4>{t)) sin 2 8(t) 

(S10) 



Mo 

where Bo e (<f)(i)) = —MgQ(N xx — N yy )sin(2(j)(t)). The torque acting on the magnetization within 
unit volume due to shape and stress anisotropy is 



T E (*) = -n m (t)xVE[8(t),cf>(t)} 

= -e r x [{2B((/>(t))sin6(t)cos9(t)}eo - {B 0e (<j)(t)) sin6{t)}e^] 

= -{2B(^(t))sin9{t)cos6{t)}e (j) - {B Oe (0(t)) sin6{t)}e e (Sll) 



The magnetization dynamics of the single-domain magnet under the action of various torques 
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is described by the Landau-Lifshitz-Gilbert (LLG) equation as follows. 



dn m (t) 



dt 



+ a\ n m (t) x 



dn m (t) \ 7 



dt 



My 



■Te(*) 



(S12) 



where a is the dimensionless phenomenological Gilbert damping constant, 7 = 2//_b/^ is the gyro- 
magnetic ratio for electrons and is given by 2.21 x 10 5 (rad.m).(A.s) -1 , and My = [i§M s Vl. In the 
spherical coordinate system, 



dn m (t) 
dt 



6'(t)e g + sin9(t)(j)'(t) e . 



(S13) 



where the prime denotes first derivative with respect to time. Accordingly, 



a 



n m (i) x 



V 



dn m (t) 
~~dT 



-asinOit) (f>'(t) e e + aO'(t) 



and 



dn m (t) 
dT 



( 



■ + a 



n m (i) x 



V 



dn m (t) 
~~dT 



(S14) 



{6'(t) - asinOit) <j>\t)) e e + (sin6{t) <j>'{t) + aO'{t)) e^. (S15) 



Equating the eg and e<h components in both sides of Equation (|S12p . we get 



9' {t) - asinO{t) <j)'(t) 
sin9{t) 4>'(t) + ad'(t) 



7 

My 

7 

My 



B 0e (<p(t)) sin6(t) 
2B(cp(t))sine(t)cos9{t). 



Simplifying the above, we get 



(S16a) 
(S16b) 



(1 + a 2 ) 6'{t) 
(1 + a 2 ) 4>'{t) 



M v 



[B Qe (4>(t))sin6{t) + 2aB(<f>(t))sin6(t)cosO(t)] 



My 



[aB 0e (<j)(t)) - 2B((p(t))cos6(t)] (sin9(t) / 0). 



(S17) 
(S18) 



We will assume that the initial orientation of the magnetization is aligned close to the -z-axis 
so that 0i n itiai = 180° — e. If e = and the magnetization is exactly along the easy axis [6 = 0°, 
180°], then no amount of stress can budge it since the effective torque exerted on the magnetization 
by stress will be exactly zero (see (jSllj) ). Such locations are called "stagnation points". Therefore, 
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we will assume that e = 1°. This is not an unreasonable assumption since thermal fluctuations can 
dislodge the magnetization from the easy axis and make e — > 1°. 

We should notice from Equation ()S17j) that there is the possibility of one more stagnation point 
at 9(t) = 4>{t) = 90° [in-plane hard axis] since there 0'(t) =0. At 9 = 90°, Equation (fS17^ becomes 

(1 + a 2 ) 9'{t) = -jL. B Qe m)) = -jL- Y M*n(N xx - N yy )sin(2^{t)) (S19) 

which indicates that as long as < 90°, the magnetization vector will continue to rotate towards 
the correct final state without being stuck at 9 = 90°. This will avoid stagnation. Note that when 9 
= 90°, we will stagnate if 4>(t) = 90°, rotate back towards the intial state along the -z-axis (wrong 
state) if 4>(t) > 90°, and rotate towards the correct state along the +z-axis if 4>(t) < 90°. 

At high enough stress, the out-of-plane excursion of the magnetization vector is significant 
and (p(t) < 90° so that stagnation is indeed avoided and the correct state is invariably reached. 
However, at low stress, the first term in Equation (|S18p will suppress out-of-plane excursion of the 
magnetization vector and try to constrain it to the nanomagnet's plane, thereby making 4>{t) = 90°. 
This will result in stagnation when 9 reaches 90° and switching will fail. Whether this happens or 
not depends on the relative strengths of the two terms in Equation (|S18p that counter each other. 
We need to avoid such low stresses to ensure successful switching. Thus, there is a minimum value of 
stress for which switching takes place. This minimum value is determined by material parameters. 

One other issue deserves mention. We have shown explicitly that we can switch from an initial 
state close to the -z-axis to a final state close to the +z-axis. Can we do the opposite and switch 
from +z-axis to -z-axis? For a single isolated magnet, this is always possible and the dynamics 
is identical. In magnetic random access memory (MRAM) systems, there are two strongly dipole 
coupled magnet in close proximity, where one is the soft magnetic layer and the other is the hard 
magnetic layer. In that scenario, there is a difference between switching from the anti-parallel to the 
parallel and from the parallel to the anti-parallel arrangement of the two magnets owing to dipole 
coupling which is different in the two cases. This is not an issue for the single isolated magnet 
considered here. 

S2 Material parameters 

The material parameters that are used in the simulation are given in the Table Q] [321~9Z! • They 
ensure that the shape anisotropy energy barrier is ~32 kT. 
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Terfenol-D 


Nickel 


Cobalt 


I\/lfl10T* Fl"V"l<i I PI 1 
IV-LdJUI Ct-A-lO Id 1 


101 / ^ n tv i 

± VJ ± . 1 tj 11111 


1(1^ n tti 

±UU 11111 


101 / ri n in 
1U1. 1 U 11111 


i n nr pivit; ( v~\\ 

1V.1111U1 CLA.10 1 VJ I 


Z/\j.AO 11111 


Q n n m 

C*(J 11111 


Q8 95 nm 

J70.aU 11111 


Thickness (t) 


10 nm 


10 nm 


10 nm 


Young's modulus (Y) 


8xl0 lu Pa 


2.14xl0 n Pa 


2.09xlO ri Pa 


Magnetostrictive coefficient ((3/2)A s ) 


+90xl0" 5 


-3xl0~ 5 


-3xl0" 5 


Saturation magnetization (M s ) 


8xl0 5 A/m 


4.84 xlO 5 A/m 


8xl0 5 A/m 


Gilbert's damping constant (a) 


0.1 


0.045 


0.01 



Table 1: Material parameters for different materials. 



S3 Procedure for determining the voltage required to generate a 
given stress in a magnetostrictive material 

In order to generate a stress a in a magnetostrictive layer, the strain in that material must be 
e = cr/Y, where Y is the Young's modulus of the material. We will assume that a voltage applied 
to the PZT layer strains it and since the PZT layer is much thicker than the magnetostrictive layer, 
all the strain generated in the PZT layer is transferred completely to the magnetostrictive layer. 
Therefore, the strain in the PZT layer must also be e. The electric field needed to generate this 
strain is calculated from the piezoelectric coefficient cfei of PZT (c^i = 1.8 x 10~ 10 m/V [S8J) and 
the corresponding voltage is found by multiplying this field with the thickness of the PZt layer. 



S4 Calculation of the energy Ed dissipated internally within the 
magnet due to Gilbert damping 

Because of Gilbert damping in the magnet, an additional energy is dissipated when the magnet 
switches. This energy is given by the expression 



o 



P d (t)dt, (S20) 
where r is the switching delay and Pd(t), the dissipated power is given by [£|9"1 SfTO] 



a 7 

f *° (HJW «»'' (S21) 

We calculate this quantity numerically and add that to the quantity CV 2 dissipated in the 
switching circuit to find the total dissipation E tota i. The results are plotted as a function of r in 
Fig. 2 of the main letter. 
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S5 Simulation results 

Some additional simulation results and corresponding discussions are given in the Figures 1521- IS51 
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Figure SI: Voltage required to switch a multiferroic nanomagnet of the shape and size considered 
in this paper versus switching delay. Three different layers are considered for the magnetostrictive 
layer. Terfenol-D requires the smallest voltage since it has the highest magnetostrictive coefficient. 
This tiny voltage requirement makes this mode of switching magnets extremely energy-efficient. 
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Dynamics of e with time (Terfenol-D, 1.92 MPa Stress) 
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Figure S2: Time evolution of the polar angle 9 in a Terfenol-D/PZT multiferroic nanomagnet un- 
der 1.92 MPa applied stress. Stress is applied abruptly at time t = to rotate the magnetization 
away from its initial orientation close to the -z-axis (9 = 179°) and it is removed abruptly once 
9 reaches 90°, which corresponds approximately to the hard axis. Thereafter, the magnetization 
spontaneously decays to the easy axis since shape anisotropy prefers the unstressed magnet's mag- 
netization to align along the easy axis. Whether it decays to the +z-axis or -z-axis is determined 
by the sign of l?oe when 9 reached 90°. If the sign is positive, then 9' is negative and 9 will decrease 
with time, finally reaching the value of 1° so that the magnetization aligns along the desired +z-axis. 
It is therefore imperative to ensure that Bo e is positive, which will happen only if <j) < 90° when 
9 = 90°. We show in the next figure that this indeed happens as a consequence of the coupled 6- 
and (^-dynamics. The coupled dynamics therefore plays a critical role to ensure correct switching. 
Note that the magnetization spends a lot of time around 9 = 90° which is the hard axis. Once the 
magnetization gets past the hard axis, it quickly reaches the easy axis. This can be understood 
by looking at the energy profile in Fig. IS5I The small stress causes a shallow energy minimum at 
the hard axis, but upon removal of stress, shape anisotropy causes a tall energy barrier at the hard 
axis. Hence it is much easier to approach the easy axis from the hard axis, but much harder to 
approach the hard axis from the easy axis. The total switching delay in this case is ~ 100 ns, out 
of which nearly ~90 ns is spent to get past the hard axis starting out from the easy axis, and ~10 
ns to decay to the easy axis from the hard axis. 
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Dynamics of * with time (Terfenol-D f 1.92 MPa Stress) 
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Figure S3: Time evolution of the azimuthal angle (ft during the switching of a Terfenol-D/PZT 
multiferroic nanomagnet subjected to 1.92 MPa stress. Initially, we start from <p = 90° and 9 = 
179° and therefore the second term in Equation (|S18p . B(<p(t))cosQ starts out as positive (since 
B{(j)(t)) goes negative upon application of stress and cos 9 is negative in the range 90° < 9 < 
180°). Consequently </> decreases with time initially. This decrease of <f> affects the Bo e ((j){t)) term 
in Equation (|S17p facilitating the rotation of the magnetization angle 9 from 179° toward 90°. 
Thereafter, (ft starts to increase because the term BQ e {(p(t)) becomes non-zero as soon as deviates 
from 90°. But cj) never reaches exactly 90° when 9 = 90°. This avoids a possible stagnation point at 
the hard axis. When 9 = 90° [hard axis], stress is removed but the finite value of Bo e ((p(t)) [the term 
due to shape anisotropy] continues to rotate the magnetization towards the +z-axis. When 9 < 90°, 
the term B{<f>(t))cos0 goes positive and according to Equation (|S18[) . <j> starts to decrease. As <f> 
decreases, Bo e ((f>(t)) increases and according to Equation (|S17p . 9 decreases sharply and ultimately 
reaches a value of 1°, at which point the switching is complete. Note that <p never deviates too far 
from 90° at this low value of stress (except at the very tail end of the switching), meaning that 
the magnetization vector is pretty much confined to the plane of the magnet (y-z plane) and its 
out-of-plane excursion is very small. Under high stresses, the out-of-plane excursion can be quite 
significant. 
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Dynamics of magnetization axis (Terfenol-D,1.92 MPa Stress) 




Figure S4: Trajectory traced out by the tip of the magnetization vector in space during switching. 
The magnet is a Terfenol-D/PZT multiferroic nanomagnet subjected to 1.92 MPa stress. 
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~ i 1 1 1 1 r 



i t 



S 



■With stress 



no stress 



/ 



/ 



\ 



\ 



\ 



\ 



v - 



J I I I I I I L 



20 40 60 80 100 120 140 160 180 

e (degree) 



Figure S5: Steady state energy profiles of a stressed and unstressed Terfenol-D/PZT multiferroic 
nanomagnet. The magnitude of the stress is 1.92 MPa. Without any applied stress, the potential 
profile depicts the shape anisotropy energy barrier which is 0.8 eV or 32 kT. 
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